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1. Introduction 

Daryl Geller and I started to work on this paper, dedicated to the memory of 
Leon Ehrenpreis, in the Fall of 2010. Sadly, Daryl Geller passed away suddenly in 
January of 2011. I will always remember him as a good friend and a wonderful 
mathematician. 

Analysis on two dimensional surfaces and in particular on the sphere S 2 found 
many applications in computerized tomography, statistics, signal analysis, seismol- 
ogy, weather prediction, and computer vision. During last years many problems of 
classical harmonic analysis were developed for functions on manifolds and especially 
for functions on spheres: splines, interpolation, approximation, different aspects of 
Fourier analysis, continuous and discrete wavelet transform, quadrature formulas. 
Our list of references is very far from being complete PQ-[B], [8]- [IB], [H]-[3S]. More 
references can be found in monographs |12j . |20j . 

The goal of the paper is to describe three types of cubature formulas on general 
compact Riemannian manifolds which require essentially optimal number of nodes. 
Cubature formulas introduced in section 3 are exact on subspaces of band-limited 
functions. Cubature formulas constructed in section 4 are exact on spaces of varia- 
tional splines and, at the same time, asymptotically exact on spaces of band-limited 
functions. In section 5 we prove existence of cubature formulas with positive weights 
which are exact on spaces of band-limited functions. 

In section 7 we prove that on homogeneous compact manifolds the product of two 
band-limited functions is also band-limited. This result makes our findings about 
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cubature formulas relevant to Fourier transform on homogeneous compact mani- 
folds and allows exact computation of Fourier coefficients of band-limited functions 
on compact homogeneous manifolds. 

It worth to note that all results of the first four sections hold true even for non- 
compact Riemannian manifolds of bounded geometry. In this case one has properly 
define spaces of bandlimited functions on non-compact manifolds |26] . 

Let M be a compact Riemannian manifold and £ is a differential elliptic op- 
erator which is self adjoint in L2(M) = L2(M, dx), where dx is the Riemannian 
measure. The spectrum of this operator, say = Ao < Ai < A2 < is discrete 
and approaches infinity. Let uq, u\, 112, ••• be a corresponding complete system of 
real- valued orthonormal eigenf unctions, and let E w (£), uj > 0, be the span of all 
eigenfunctions of C, whose corresponding eigenvalues are not greater than u). For 
a function / £ L2CM.) its Fourier transform is the set of coefficients {cj-(/)}, which 
are given by formulas 



Jm 

By a discrete Fourier transform we understand a discretization of the above formula. 
Our goal in this paper is to develop cubature formulas of the form 



where {xk} is a discrete set of points on M and {wk} is a set of weights. When 
creating such formulas one has to address (among others) the following problems: 

1. to make sure that there exists a relatively large class of functions on which 
such formulas are exact; 

2. to be able to estimate accuracy of such formulas for general functions; 

3. to describe optimal sets of points {xk} for which the cubature formulas exist; 

4. to provide "constructive" ways for determining optimal sets of points {x^}; 

5. to provide "constructive" ways of determining weights Wk', 

6. to describe properties of appropriate weights. 

In the first five sections of the paper we construct cubature formulas on general 
compact Riemannian manifolds and general elliptic second order differential oper- 
ators. Namely, we have two types of cubature formulas: formulas which are exact 
on spaces E w (£) (see section 3), i. e. 



and formulas which are exact on spaces of variational splines (see section 4) . More- 
over, the cubature formulas in section 4 are also asymptotically exact on the spaces 
E u (£). For both types of formulas we address first five issues from the list above. 
However, in the first four sections we don't discuss the issue 6 from the same list. 

In section 5 we construct another set of cubature formulas which are exact on 
spaces E w (£) which have positive weights of the "right" size. Unfortunately, for 
this set of cubatures we unable to provide constructive ways of determining weights 

If one considers integrals of the form (|I.Ij) then in the general case we do not have 
any criterion to determine whether the product fuj belongs to the space E w (£) in 



(1.1) 




(1.2) 




(1.3) 
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order to have an exact relation 



(1.4) / fuj = y2f(xk)uj(x k )w k 



for cubature rules described in sections 1-4. However, if M is a compact homoge- 
neous manifolds i. e. M = G/K, where G is a compact Lie group and K is its 
closed subgroup and C is the second order Casimir operator (see (|6.2p below) then 
we can show that for f,g£ E w (£) their product fg is in Eid u (C), where d — dimG 
(see section 7). 

2. PLANCHEREL-POLYA-TYPE INEQUALITIES 

Let B(x, r) be a metric ball on M whose center is x and radius is r. The following 
important lemma can be found in [26) . [29] . 

Lemma 2.1. There exists a natural number Njsa, such that for any sufficiently 
small p > 0, there exists a set of points {y u } such that: 

(1) the balls B(y u ,p/A) are disjoint, 

(2) the balls B{y v ,p/2) form a cover o/M, 

(3) the multiplicity of the cover by balls B(y Vl p) is not greater than JVm- 

Definition 1. Any set of points M p = {y„} which is as described in Lemma [2.11 
will be called a metric p-lattice. 

To define Sobolev spaces, we fix a cover B = { £?(?/„, ro)} of M of finite multi- 
plicity TVm (see Lemma |2~T|) 

(2.1) M = {J B(y v ,r ), 

where B{y v ,ro) is a ball centered at y v £ M of radius ro < pm, contained in a 
coordinate chart, and consider a fixed partition of unity ^ = {ip u } subordinate to 
this cover. The Sobolev spaces H S (M.), s £ R, are introduced as the completion of 
C°°(M) with respect to the norm 

1/2 

(2.2) 



I/IIh*(m)= ^Uuf\\%. 



(B(y u ,r )) 



Any two such norms are equivalent. Note that spaces H S (M.), s £ R, are domains 
of operators A s / 2 for all elliptic differential operators A of order 2. It implies, that 
for any s £ R there exist positive constants a(s), b(s) (which depend on A) such 
that 

1 /2 

(2.3) ||/||^ ( M) < a(s) (||/||i 2(M) + P s/2 /Hl 2( m)) < b(s)\\f\\ H , m 
for all / £ H S (M). 

We are going to keep notations from the introduction. Since the operator C is of 
order two, the dimension Af u of the space E w (£) is given asymptotically by Weyl's 
formula 

(2.4) A/L(M) x C(M)w™/ 2 , 
where n = dirnM.. 

The next two theorems were proved in [26j , [30] , for a Laplace-Beltrami operator 
in Z/2(M) on a Riemannian manifold M of bounded geometry, but their proofs 
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go through for any elliptic second-order differential operator in Z/2(M). In what 
follows the notation n — dim M is used. 

Theorem 2.2. There exist constants C\ > and po > 0, such that for any natural 
m > n/2, any < p < po, and any p-lattice M p = {xk} 7 the following inequality 
holds: 

1/2 

(2-5) ( V |/(^)| 2 ) <C lp -"/2||/|| ffm(M)) 



y x k eM p 



for all f E H m (M). 



Theorem 2.3. There exist constants Ci > 0, and po > 0, such that for any natural 
m > n/2, any < p < po, and any p-lattice M p = {xk} the following inequality 
holds 

(2.6) ||/||h«(m) < C 2 { P n ' 2 ( E l/(^)l 2 ) +P m \\£ m/2 f\\L 2m 

As one can easily verify the norm of C on the subspace E w (£) (the span of 
eigenfunctions whose eigenvalues < lj) is exactly uj. In particular one has the 
following Bernstein-type inequality 

(2-7) l|£'/IU a ( M ) <^II/IIl 2( m), seR+, 

for all / € E U (C). This fact and the previous two theorems imply the following 
Plancherel-Polya-type inequalities. Such inequalities are also known as Marcinkewicz- 
Zygmund inequalities. 

Theorem 2.4. Set mo = [S-l +1. If C\, Ci are the same as above, a(mo) is from 

\2. 3\) , and cq = (IC^ 1 ) 1 '"" 1 then for any tu > 0, and for every metric p-lattice 
M p = {x^ with p = cqlu~ x / 2 , the following Plancherel-Polya inequalities hold: 

C 1 - 1 fl(m )- 1 (l + ^)- rao/2 (El/^)l 2 ) <P~ n/a ll/lk(M) < 

(2-8) (2C 2 )(^2\f(x k )\^ , 

for all f € E w (£) and n = dim M. 

Proof. Since C is an elliptic second-order differential operator on a compact man- 
ifold which is self-adjoint and positive definite in L 2 (M) the norm on the Sobolev 
space iJ mo (M) is equivalent to the norm ||/||l 2 (m) + l|£ mo/2 /llL 2 (M)- Thus, the 
inequality (|2.5I) implies 

v 1/2 

E l/(^)| 2 <Ci a (moV-" /2 (||/||L 2 (M) + ||/: mo/2 /l|L 2 (M)). 

,x k eM p J 
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The Bernstein inequality shows that for all / € E u (£) and all lu > 

H/llL 2 (M) + l|/: mo / 2 /llL 2 (M) < (l + ^r n/2 ||/|U 2 (M). 

Thus we proved the inequality 

(2-9) ^ 

Cr 1 a(m )- 1 (l+ W )-"V2 | £ \f(x k )\A <p-"^\\f\\ L2m , /€E U (£). 

To prove the opposite inequality we start with inequality (|2.6p where m = 
[§] + 1. Applying the Bernstein inequality (|2.7[) we obtain 

(2-10) ||/|U 2 (M)<CV l/2 ( E l/(^)| 2 ) +C 2 p m ^ mo/2 ||/||L 2 (M), 

where / 6 E u (£). Now we fix the following value for p 

/i \ l/mo /, \ l/mo 

With such p the factor in the front of the last term in (12.10)) is exactly 1/2. Thus, 
this term can be moved to the left side of the formula (12.10)) to obtain 

(2.H) i|i/iu 2( M)<c 2 p"/ 2 ( \n*k)A ■ 

\x k eM p J 

In other words, we obtain the inequality 

P-" /2 ||/IIl 2 (m)<2C 2 I Yl l/(**)| 2 ' 

\x k £M p 

The theorem is proved. 

□ 

It is interesting to note that our p-lattices (appearing in the previous Theorems) 
always produce sampling sets with essentially optimal number of sampling points. 
In other words, the number of points in a sampling set for E w (£) is "almost" the 
same as the dimension of the space E w (£) which is given by the Weyl's formula 
(El. 



Theorem 2.5. // the constant Co > is the same as above, then for any lu > 
and p = cqlo^ 1 / 2 , there exist positive at, a 2 such that the number of points in any 
p-lattice M p satisfies the following inequalities 

(2.12) a X Lu n ' 2 < \M p \ < a 2 Lo n/2 ; 

Proof. According to the definition of a lattice M p we have 

\M p \ inf Vol(B{x,p/4)) < Vol(M) < \M p \ sup Vol(B{x, p/2)) 

xeM xeM 

or 

« < |MJ < ™M . 

sup, e M Vol(B(x, p/2)) " 1 pl " inf^M Vol (B{x, p/A)) 
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Since for certain ci(M),C2(M), all x € M and all sufficiently small p > one 
has a double inequality 

ci(M)p" < Vol(B(x,p)) < c 2 {M)p n , 

and since p = c cj -1 / 2 , we obtain the inequalities (|2.12[) for certain a\ — a^M), a 2 = 
a 2 (M). 

□ 

3. CUBATURE FORMULAS ON MANIFOLDS WHICH ARE EXACT ON BAND-LIMITED 

FUNCTIONS 



Theorem 12.41 shows that if Xk is in a p lattice M p and -dk is the orthogo- 
nal projection of the Dirac measure 5 Xk on the space E w (£) (in a Hilbert space 
ff~ n / 2 ~ e (M), e > 0) then there exist constants ci = Ci(M,£,u>) > 0, c 2 = 
c 2 (M, £) > 0, such that the following frame inequality holds for all / € E w (£) 

(3.i) d (ek.wi 2 ) <p- n/a n/iu a (M) <c 2 (ek/^>i 2 ) . 

where 

{f,0 k ) = fan), /eW 

From here by using the classical ideas of Duffin and Schacffcr about dual frames 
[7] we obtain the following reconstruction formula. 

Theorem 3.1. If M p is a p-lattice in Theorem \2.4\ with p = c^lo^ 1 / 2 then there ex- 
ists a frame {Qj} in the space E u (£) such that the following reconstruction formula 
holds for all functions in E w (£) 

(3-2) /= E f(*k)®k. 

x k eM„ 



This formula implies that for any linear functional F on the space E w (£) one 
has 

F(f)= E fMF(e k ), /GE U (£). 

x k GM p 

In particular, we have the following exact cubature formula. 

Theorem 3.2. If M p is a p-lattice in Theorem \2.J\ with p = cqoj^ 1 ' 2 and 

v k = f Ok, 

JTsA 

then for all f £ E u (£) the following holds 

(3-3) f f= E /foK. /GE-OC). 

Thus, we have a cubature formula which is exact on the space E w (£). Now,we 
are going to consider general functions / e L 2 (M). Let f u be orthogonal projection 
of / onto space E w (£). As it was shown in [52] there exts a constant Ck, m that the 
following estimate holds for all / € L 2 (M) 

(3-4) \\f-f u \\L a <M<^n m -k(£ k f,l/u), fc,meN. 
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Here the modulus of continuity is denned as 
(3.5) 



n r (g,s)= sup ||A; 5 ||, g£L 2 (M), r £ N, 

\t\<s 



where 
(3.6) 



r 

A r T g = (-iy +1 J2(-^ 1 C ] r e^g, r £ R, r£ 

3=0 



N. 



Thus, by combining (|3.3j) and (|3.4|) we obtain the following theorem. 

Theorem 3.3. There exists a Cq = Co(M, C) and for any < k < to, k, m £ N, 
there exists a constant Ck, m > such that if M p = {xk} is a p-lattice with < p < 
cqlo^ 1 then for the same weights {vj} as in 



(3.7) 



where f u is the orthogonal projection of f £ L 2 (M) onto E w (£). 

Note (see [32]), that / £ L2(M) belongs to the Besov space B^^M) if and 
only if 

Q m (f,l/u;) = 0(u J - a ), 

when uj — > oo. Thus, we obtain that for functions in B 2 oo(M) the following 
relation holds 



(3.8) 



/ /-yWw(>j) 

JM 



0{uj- a ), uj — >oo. 



4. CUBATURE FORMULAS ON COMPACT MANIFOLDS WHICH ARE EXACT ON 

VARIATIONAL SPLINES 



Given a p lattice M p = {x 7 } and a sequence {z 7 } £ l 2 we will be interested to 
find a function Sk £ H 2k (M.), where k is large enough, such that 

(1) Sfe(x 7 ) = z 1 ,x 1 £ M p ; 

(2) function Sk minimizes functional g — > ||£ fe g||L 2 (M)- 

We already know (|2.5[) , (|2.6p that for k > d the norm H 2k (M) is equivalent to the 
norm 

Cx{p)\\f\\ H * H M) < \\c k f\\ L2{ M) + £ l/^)! 2 < c 2 (p)\\f\\ H 2k {M.)- 

\x^£M p J 

For the given sequence {z 7 } £ l 2 consider a function / from H 2k (M) such that 
f{x-y) = Zj. Let Pf denote the orthogonal projection of this function / in the 
Hilbert space H 2k (M.) with the inner product 

</,<?>= E /WsK)+<^ 2 /,^ /2 !?> 
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on the subspace U 2k (M p ) = {/ G H 2k (M)\f(x 1 ) — 0} with the norm generated by 
the same inner product. Then the function g = / — Pf will be the unique solution 
of the above minimization problem for the functional g ||£ fe g||L 2 (M)j k > d. 
Different parts of the following theorem can be found in [5Tj . 

Theorem 4.1. The following statements hold: 

(1) for any function f from H 2k (M.), k > d, there exists a unique function 
Sk(f) from the Sobolev space H 2k (M.), such that /|m p = Sfe (/) |jvf„ ; and this 
function Sk(f) minimizes the functional u — > \\C k u\\^ 2 (^jiy, 

(2) every such function Sk(f) is of the form 

s k (f)= ^ <•'■• 



where the function L 2 G H 2 (M), x-y G M p minimizes the same functional 
and takes value 1 at the point x 7 and at all other points of M p ; 

(3) functions L 2k form a Riesz basis in the space of all polyharmonic functions 
with singularities on M p i.e. in the space of such functions from H 2k (M.) 
which in the sense of distributions satisfy equation 

C 2k u= a 1 5{x 1 ) 

where 5(Xy) is the Dirac measure at the point x 7 ; 

(4) if in addition the set M p is invariant under some subgroup of diffeomor- 
phisms acting on M then every two functions L 2k , L 2k are translates of 
each other. 

The crucial role in the proof of the above Theorem 14.11 belongs to the following 
lemma which was proved in |26j . 

Lemma 4.2. A function f 6 La(M) satisfies equation 

C 2k f= £ a^( Xj ), 

where {a 7 } € h if and only if f is a solution to the minimization problem stated 
above. 

Next, if / £ H 2k (M), k > d, then / - s k (f) G U 2k (M p ) and we have for k > d, 

11/ - ^(/)IU 2 (M) < (C 0P ) k \\£ k/2 (f - fl fc (/))||i a (M). 

Using minimization property of Sfc(/) we obtain the inequality 



(4.1) 



/- XI f( x i) L * 



< ( C0 p) fe ||£ fe/2 /ll 



L 2 (M), 



k > d. 



L 2 (M) 



and for / G E w (£) the Bernstein inequality gives for any / G E w (£) 



(4.2) 



< ( Co pv^) fc ||/ll 



L 2 (M) 



L 2 (M) 



for k > d. The last inequality shows in particular, that for any / G E w (£) one has 
the following reconstruction algorithm. 
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Theorem 4.3. There exists a cq = co(M) such that for any u> > and any M p 

with p = cqlu^ 1 the following reconstruction formula holds in L2(M)-norm 

(4.3) /= lim V f( Xj )L<V, k>d, 

l— >oo — * J 

for allfe-E u (C). 

To develop a cubature formula we introduce the notation 



(4.4) 



^7 



I L^(x)dx, 

M 



where L x G S k (M p ) is the Lagrangian spline at the node x 1 . 
Theorem 4.4. (1) For any f e H 2k (M) one has 



(4.5) 



(4.6) 



/ fdxn J2 A ? } /(^), k>d, 



and the error given by the inequality 



f fdx- A f/(^7) 



< ^(M)( Co p) fc ||£ fc / 2 /|| 



L 2 (M) 



for k > d. For a fixed function f the right-hand side of 114- 6\ ) goes to zero 
as long as p goes to zero. 
(2) The formula j4-5\ ) is exact for any variational spline f g S k (M p ) of order 
k with singularities on M p . 

By applying the Bernstein inequality we obtain the following theorem. This 
result explains our term " asymptotically correct cubature formulas" . 

Theorem 4.5. For any f g E w (£) one has 



(4.7) 



/ fdx- J2 A 7 fe) /(z 7 ) 



< Vol(M)(c pV^) k \\f\\ 



L 2 (M) 



for k > d. If p = cquj 1 / 2 the right-hand side in \4- 7p goes to zero for all f G E w (£) 
as long as k goes to infinity. 

5. Positive cubature formulas on compact manifolds 

Let M p = {xk}, k = 1, N(M p ), be a p-lattice on M. We construct the Voronoi 
partition of M associated to the set M p = {xk}, k = 1, N(M p ). Elements of this 
partition will be denoted as Mk,p- Let us recall that the distance from each point 
in Mj tP to Xj is less than or equal to its distance to any other point of the family 
M p = {xk} 1 k — 1, ...,N(M p ). Some properties of this cover of M are summarized 
in the following Lemma, which follows easily from the definitions. 



1. 



,N(M p ), have the following properties: 



Lemma 5.1. The sets Aik, P , k 

1 ) they are measurable; 

2) they are disjoint; 

3) they form a cover of M; 

4) there exist positive a\, a^, independent of p and the lattice M p — {xk}, such 
that 
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(5.1) 



aip" < fi(Mk, P ) < a 2 p n - 



In what follows we are using partition of unity VP = {ipv} which appears in (|2.2|) . 
Our next goal is to prove the following fact. 

Theorem 5.2. Say p > 0, and let {M k ,p\ be the disjoint cover of M which is 
associated with a p-lattice M p . If p is sufficiently small then for any sufficiently 
large K G N there exists a C{K) > such that for all smooth functions f the 
following inequality holds: 



Y Y ^f( x k) VM k , p - / f{x)d 



< 



(5.2) 



v x k £M p 

C(K)f2p n ^\\(I + ^ m/2 f\U M ), 

l/3|=l 

where C(K) is independent of p and the p-lattice M p . 
Proof. We start with the Taylor series 

(5.3) i> v f{y)-^ v f{x k )= V ±d a (^f)(x k )(x k - y r+ 

\<\a\<m-\ 

j2 - r t m - x d a i> v f{x k + td)6 a dt, 

\a\—m 

where / G C°°(R d ), y G B(x k ,p/2), x = (x^\...,x^), y = (yW, yW), a = 
( ai ,...,a d ), (x-y) a = (xW-yW) ai ...(xW-yW) a *, r = \\x- Xi % 9 = (x-x^/t. 

We are going to use the following inequality, which is essentially the Sobolev 
imbedding theorem: 

(5.4) \U>vf)(x k )\ < C n , m y ^ n/P \\^f)\\w l( B( Xk , P ) V 1 ^ P ^ °°. 

0<j<m 

where m > n/p, and the functions {?/v} form the partition of unity which we used 
to define the Sobolev norm in (|2.2j) . Using (|5.4j) for p = 1 we obtain that the 
following inequality 



< 



(5.5) Y -,d a {^f)(x k ){x k - V y 

\<\a\<m-\ 

C(n,m)pM Y Y P b ^ n \\d a+1 (^f)\\L l{ B {xk , p)) , m>n, 

for some C(n, m) > 0. Since, by the Schwarz inequality, 

(5.6) ||9 q (Vv/)I|l i( s(^p)) < C(n)p n / 2 \\d a (^f)\\ L2{B{xk , p)) 
we obtain the following estimate, which holds for small p: 



(5.7) 



sup 

y£B(x k ,p) 



Y -,d a {^f){x k ){x k - y y 



l<\a\<m-l 



< 
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C{n,m) P m ~ n/2 W dP (^f)h 2 {B( Xk , p) ), m>n. 

l<|/3|<2m 

Next, using the Schwarz inequality and the assumption that m > n — dim M, \a\ — 
to, we obtain 

t m - 1 d a ^f( Xk + te)0 a dt < 



1/2 



t m - n / 2 - 1 / 2 \t n l 2 - 1 l 2 d a ij) l/ f{x k +t6)\dt < 



C t 



4.2m— n— 1 



1/2 



t n - x \d a ^ v f{x k + t6)\ 2 dt) ■ 

1/2 



m > n. 



We square this inequality, and integrate both sides of it over the ball B(x k , p/2), 
using the spherical coordinate system (t, 0) . We find 



C(m, n) 



B{x k ,p) 
p/2 





p/2 

C(m,n) I t 



t m - l d a Tp v f{x k + t9)9 a dt 

271 



T n ~ dOdT < 







t n - 1 d a (^f){ Xk + te)9 a dt 







T n - Y dQdr < 



27T /./j/2 



2m — n I aa 



|a Q (i/' l //)(a; fe +t6>)| z r n - 1 dTd6»j eft < 

C m , nP 2 l Q l||9«(^/)|| 2 L2(B(;Cfc , p)) , 

where t = [| ac — Xk\\ < p/2, m = \a\ > n. Let {Mk,p} be the Voronoi cover of M 
which is associated with a p-lattice M p (see Lemma ISTTj) . From here we obtain 



(5.8) 



H>uf{y) - 4>vf{x k )\dx < 



C(n,m) X! ^ l+n/a H^(^/)lk(B 



K|fl|<2r, 



y 1 



t m - 1 d a ^f(x k + t9)6 a dt 



< 



\a\—m \ v 



C(n,m) J] P m+n/ V(M)hz(B(* k , P ))+ 

l<|/3|<2m 

2 \ Va 

T n ~ X dTdd 



p) 



t m - i d a ^f{x k + te)9 a dt 



C(n,m) J2 P M+n/ V(^f)h 2mXh , P )). 

l<\0\<2m 

Next, we have the following inequalities 



Y ^f( x k) P-M k , P - / f{x)d. 
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H M Tp y f(x)dx ~ ^ ^f( x k) (J>M k>p < 
v \ k J Mh, P k J 



(5.9) 



EE 

v k 



ip„f{x) - ip u f{x k ) pM k , P dx 



<C(n,m)p"/ 2 £ E E P m \\d^^f)\\L 2 (B( Xk 



p))> 



v x h EM„ K|fl|<2r, 



where m > n. Using the definition of the Sobolev norm and elliptic regularity 
of the operator I + C, where I is the identity operator on Z,2(M), we obtain the 
inequality (|5.2p . □ 



Now we are going to prove existence of cubature formulas which are exact on 
E W (M), and have positive coefficients of the "right" size. 

Theorem 5.3. There exists a positive constant ao, such that if p = ao(uj + l) -1 / 2 , 
then for any p-lattice M p = {xk}, there exist strictly positive coefficients p Xk > 
0, Xk € M p , for which the following equality holds for all functions in E w (£): 



(5.10) 



/ fdx= Y Vx k f(xk)- 



Moreover, there exists constants a, c 2 , such that the following inequalities hold: 
(5.11) cip n < p Xk < c 2 p n , n = dim M. 

Proof. By using the Bernstein inequality, and our Plancherel-Polya inequalities 
(|2.4[) . and assuming that 



(5.12) 



P < 



1 



2VW+T 

we obtain from (|5.2p the following inequality: 



E E ^f( x k)pM k , P - f{x)d. 

v x k £M, Jm 



A' 



J/31 



f\\L 2 ( 



M 



< 



1/2 



(5.13) 



C 2P n ( P VT+^ £ \f(xkt 



yX k £M p 



where C 2 is independent of p e (0, (2y/uj + l) 1 ) and the p-latticc M p . 

Let Ruj(£) denote the space of real- valued functions in E w (£). Since the eigen- 
functions of C may be taken to be real, we have E w (£) = R u {£) + iRu(C), so it is 
enough to show that (|5.10l) holds for all / € R U (L). 

Consider the sampling operator 

S:f-> {f{x k )}x k eM P , 

which maps i?^(£) into the space Rl^pl with the £ 2 norm. Let V — S(R L0 (C)) be 
the image of Ruj(£) under S. V is a subspace of M} Mp ^, and we consider it with the 



CUBATURE FORMULAS AND DISCRETE FOURIER TRANSFORM ON MANIFOLDS 13 



induced i 2 norm. If u £ V, denote the linear functional y — > (y,w) on V by i u . By 
our Plancherel-Polya inequalities (|2.4[) , the map 



{f{xk)}x h eM. 



fdx 



M 



is a well-defined linear functional on the finite dimensional space V, and so equals 
£ v for some v € V , which may or may not have all components positive. On the 
other hand, if w is the vector with components {v>(-M-k, P )} , %k € M p , then w might 
not be in V , but it has all components positive and of the right size 

a x p n < p(M k , P ) < a 2 p n , 

for some positive a±, a 2 , independent of p and the lattice M p = {xj~\. Since, for 
any vector u € V the norm of u is exactly the norm of the corresponding functional 
£ u , inequality ([5]) tells us that 



(5.14) 



\\Pw - v\\ < \\w - v\\ < C 2 p n (pVTTui) 



where P is the orthogonal projection onto V. Accordingly, if z is the real vector 
v — Pw, then 

(5.15) v + (I - P)w = w + z, 

where ||z| < C 2 p n (pV^ + w) ■ Note, that all components of the vector w are 
of order 0(p n ), while the order of ||z|| is 0(p n+1 ). Accordingly, if py/l + uj is 
sufficiently small, then p := w + z has all components positive and of the right size. 
Since p = v + (I — P)w, the linear functional y —> (y,p) on V equals l v . In other 
words, if the vector p has components {p Xk }, %k & M p , then 



E /( 



Xk)Px 



x k EM„ 



fdx 



M 



for all / € i?^(£), and hence for all / G E w (£), as desired. 



□ 



We obviously have the following result. 

Theorem 5.4. (1) There exists a cq = cq(M, C) and for any < k < m, k, m £ 

N, there exists a constant Ck, m > such that if M p — {xk} is a p-lattice 
with < p < cquj^ 1 then for the same weights {p Xj } cts in \5.10\l 



(5.16) 



/ f -y2L(xj)p X] 



(2) For functions in B2 oo(M) the following relation holds 



(5.17) 



0(cu- a ), uj — >oo. 



where f u is the orthogonal projection of f G L 2 {M.) onto E w (£). 
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6. Harmonic Analysis on Compact homogeneous manifolds 

We review some very basic notions of harmonic analysis on compact homoge- 
neous manifolds [T7], Ch. II. 

Let M, dim M = n, be a compact connected C7°°-manifold. One says that a 
compact Lie group G effectively acts on M as a group of diffcomorphisms if: 

1) every element g G G can be identified with a diffeomorphism 

g : M -> M 

of M onto itself and 

gig 2 ■ x = gi ■ (g 2 ■ x), gi,g 2 eG, x G M, 

where gig 2 is the product in G and g ■ x is the image of x under g, 

2) the identity e £ G corresponds to the trivial diffeomorphism 

(6.1) e ■ x — x, 

3) for every g G G, g ^ e, there exists a point x G M such that g ■ x =/= x. 

A group G acts on M transitively if in addition to 1)- 3) the following property 
holds: 

4) for any two points ijgM there exists a diffeomorphism g £ G such that 

9 ■ x = y. 

A homogeneous compact manifold M is a C°°-compact manifold on which a 
compact Lie group G acts transitively. In this case M is necessarily of the form 
G/K, where K is a closed subgroup of G. The notation /^(M), is used for the 
usual Banach spaces L 2 (M., dx), where dx is an invariant measure. 

Every element X of the (real) Lie algebra of G generates a vector field on M, 
which we will denote by the same letter X. Namely, for a smooth function / on M 
one has 

w t^o t 

for every x G M. In the future we will consider on M only such vector fields. The 
translations along integral curves of such vector fields X on M can be identified 
with a one-parameter group of diffeomorphisms of M, which is usually denoted as 
exptX, — oo < t < oo. At the same time, the one-parameter group exptX, — oo < 
t < oo, can be treated as a strongly continuous one-parameter group of operators 
acting on the space L 2 (M.). These operators act on functions according to the 
formula 

/ -> f(exptX -x), isR, / e L 2 (M), x G M. 

The generator of this one-parameter group will be denoted by Dx , and the group 
itself will be denoted by 

e tDx f(x) = /(exptX -x), t G M, / G L 2 (M), xeM. 

According to the general theory of one-parameter groups in Banach spaces, the 
operator Dx is a closed operator on every L 2 (M). 

If g is the Lie algebra of a compact Lie group G then ([L7], Ch. II,) it is a direct 
sum g = a+ [g, g], where a is the center of g, and [g, g] is a semi-simple algebra. 
Let Q be a positive-definite quadratic form on g which, on [g, g] , is opposite to the 
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Killing form. Let X\, ...,Xd be a basis of g, which is orthonormal with respect to 
Q. Since the form Q is yhi(G)-hivariant, the operator 

-X\ - X\ - ... - Xj, d = dim G 

is a bi-invariant operator on G. This implies in particular that the corresponding 
operator on i2(M) 

(6.2) C = -D\ -D\- ... - D 2 dl Dj = D Xj , d = dim G, 

commutes with all operators Dj = Dx r This operator C, which is usually called 
the Laplace operator, is elliptic, and is involved in most of the constructions and 
results of our paper. 

In the rest of the paper, the notation ID) = {D\, Dd}, d = dim G, will be used 
for the differential operators on X2(M), which are involved in the formula (|6.2j) . 

There are situations in which the operator C is, or is proportional to, the Laplace- 
Beltrami operator of an invariant metric on M. This happens for example, if M is 
a n-dimensional torus, a compact semi-simple Lie group, or a compact symmetric 
space of rank one. 

7. ON THE PRODUCT OF EIGENFUNCTIONS OF THE CASIMIR OPERATOR C ON 
COMPACT HOMOGENEOUS MANIFOLDS 

In this section, we will use the assumption that M is a compact homogeneous 
manifold, and that C is the operator of (I6.2p . in an essential way. 

Theorem 7.1. //M = G/K is a compact homogeneous manifold and C is defined 
as in i6.2\) , then for any f and g belonging to E w (jC), their product fg belongs to 
~&4dui(£), where d is the dimension of the group G. 

Proof. First, we show that if for an / 6 ^(M) and a positive uj there exists a 
constant C(f, lu) such that the following inequalities hold 

(7-1) l|£ fc /IU 2 (M) <C(/, W V fc ||/|| i2(M ) 

for all natural k then / € E tlJ (£). Indeed, assume that 

Xn < lu < A m +i 

and 

oo 

(7-2) / = 

°j(f) =< f> u 3 >= / f(x)uj(x)dx. 
Then by the Plancherel Theorem 

oo oo 

a^ +1 £ m 2 < £ |A^| 2 <||/: fc /lli 2 (M)<^ 2fe ll/llL(M), c = c(f, u ), 

which implies 

oo / \ 2k 

j—m-\-l x 7 
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In the last inequality the fraction u;/A m +i is strictly less than 1 and k can be 
any natural number. This shows that the series (|7.2|) docs not contain terms with 
j > m + 1, i.e. the function / belongs to E w (£). 

Now, since every smooth vector field on M is a differentiation of the algebra 
C°°(M), one has that for every operator Dj, 1 < j < d, the following equality holds 
for any two smooth functions / and g on M: 

(7-3) D j (fg) = fD j g + gD j f, 1 < j < d. 

Using formula (|6 . 2[) one can easily verify that for any natural k € N, the term 
£ k if 9) is a sum °f d k , (d = dimG), terms of the following form: 

(7.4) (/<?), i • ./: h • '/• 

For every Dj one has 

D%fg) = f(D]g) + 2(D j f)(D j g) + g(D*f). 
Thus, the function C k (fg) is a sum of (4d) k terms of the form 

(D H ...D lm f)(D 3l ...D j2k _ m g). 

This implies that 

(7.5) \C k (fg)\ < (4d) k sup sup \D h ...D im f(x)\ \D h ...D hk _ m g(y)\ . 

0<m<2fc i,jEM 

Let us show that the following inequalities hold: 

(7-6) \\D tl ...D tm f\\ L2m < ^/ 2 ||/|U 2(M ) 

and 

(7-7) \\D 3l ...D 32k _ m g\\ L2{M) < ^ 2k - m ^ 2 \\g\\ L2im 

for all f,g E E w (£). First, we note that the operator 

-£ = Df + ... + D 2 d 

commutes with every Dj (see the explanation before the formula (|6.2p ) . The same 
is true for C 1 / 2 . But then 

l|£ 1/2 /lli 2( M) =< £ 1/2 /- £ 1/2 / >=< cf, f >= 

d d d 

-£ < / >= E < ^/ >= E Ii^/IIL(M)> 
i=i j =i j= i 

and also 

l|£/||i 2(M) = I|£ 1/2 /: 1/2 /I|| 2( m) = E II^/VIILcm) = 

3=1 

Eii £l/2 A-/iiL(M)= E Pi^/IIL(M)- 
j=i j,fe=i 

From here by induction on s £ N one can obtain the following equality: 
(7.8) l|£ s/2 /llL ( M) = E llAi— A»/||i 2 (M)! s G N, 

l<ii , . . . ,i s <d 
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which implies the estimates (|7.6j) and (|7.7|) . For example, to get (|7.6j) we take a 
function / from E w (£), an m 6 N and do the following 

lA.-A^lkcM) < I E IIAv^/IILcm)' 

yi<u j m <rf 

(7.9) ll>c m/2 /llL 2 (M)<^ m/2 ll/IU 2 (M). 



dim M 
s > • 



In a similar way we obtain (|7.7|) . 

The formula (|7.5p along with the formula ([7]) imply the estimate 

l|£ fe (/<?)lk(M) < (4rf) fe sup ||A 1 ...A m /||L 2 (M)llA 1 -A 2t - m 5lU < 

0<m<2/s 

(7.10) (4rf) fe c m / 2 ||/|| L2(M ) sup HAi-A^slloo. 

0<m<2fc 

Using the Sobolev embedding Theorem and elliptic regularity of C, we obtain for 
every s > 

HAi-A^-^IU < C(M)||A 1 -A 2fc - m 9lU=(M) < 

(7.11) c(M) {llAx-A^-^slU.CM) + II^Av-A.^sIU^m)} , 

where _ff s (M) is the Sobolev space of s-regular functions on M. Since the operator 
C commutes with each of the operators A/, the estimate (|7|) gives the following 
inequality: 

\\Dh-D hk . m g\\oo < C(M){oj k - m / 2 \\g\\ L2{M) +c fc - m / 2+s || 5 |U 2(M) } < 

(7.12) 

C(M)o; fc —/ 2 {|| 5 |U 2{M) +o; s / 2 || ff || i2(M) } = C(M, a;, S )^ fc - m / 2 , 
Finally we have the following estimate: 

CI 7 Ti l IN/I 

(7.13) \\C k ifg)\\ L2m <C(MJ,g,Lu, S )(4dLu) k , s > , e N, 

which leads to our result. The Theorem is proved. □ 
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